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We calculate the singlet-triplet splitting for a couple of two-dimensional electrons in the potential 
of two positively charged impurities which are located out of plane. We consider different relations 
between vertical distances of impurities hi and /i2 and their lateral distance R. Such a system 
has never been studied in atomic physics but the methods, worked out for regular two-atomic 
molecules and helium atom, have been found to be useful. Analytical expressions for several different 
limiting configurations of impurities are obtained an interpolated formula for intermediate range 
of parameters is proposed. The i?-dependence of the splitting is shown to become weaker with 
increasing fti, /12. 

PACS: 71.70.-d, 75.30.E, 71.10.-w 
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Two-particle exchange interaction is the main mechanism of the spin-spin interaction in the impurity band of 
semiconductors. The recent experiments KM on the metallic phase in two-dimensional system show that a magnetic 
field parallel to the plane of the two-dimensional electron gas can destroy this phase. This suggests that the spin 
interaction plays an important role in the formation of the metallic state. The study of low-dimensional spin systems 
is usually based on the Heisenberg Hamiltonian J2(ik) Jik^i ' ^k7 where s is the spin 1/2 operator, and i,k denote 
different spins. Thus, the knowledge of the exchange constants Jik, which are equal to half of the singlet-triplet 
energy splitting for two electrons located at different sites is very important. 

In our previous paper we have calculated these parameters and have obtained analytical expression for the 
ground state and excited states of the Heisenberg Hamiltonian at small impurity density. We assumed, however that 
impurities are located in the same plane as electron gas. This is not the case in usual semiconductor structures, where 
impurities are separated from electrons by the so-called spacer layer and residual impurities are also out of plane 
containing electron gas. In this paper we concentrate on the problem of spin splitting in the system of two electrons 
bound to the Coulomb centers which are located outside the plain. 

Fig. 1. illustrates the general geometry of the problem. 
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FIG. 1. General geometry of the problem 



The total Hamiltonian of our artificial molecule is 



H 



Ai 
2 



A2 
2 



2 

E 



/i2 + (.X, - {-lyaf + 2;2 



1 

?'12 



^R2 + (/J2 _ /,,,)2 



(1) 



Hereafter we use effective atomic units (a.u.) which means that all distances are measured in units of the effective 
Bohr radius as = fi e/m*e'^, and energies in units of m*e^ /fi e^, where m* is the effective mass of the electron, and e 
is the dielectric constant. The projection of the distance between impurities onto the plane of electron gas is i? = 2a 
and the x-coordinates of the impurities are ±a. 

Three parameters, which determine the behavior of the system, are hi, /12, and R. Of course, there is no analytical 
solution for J {hi, /i2, R) for the whole range of parameters. In the present work we investigate the different analytical 
limits of the Hamiltonian (|^) and present an interpolated formula, which connects these limits. 

If /ii ~ /i2 ~ ^ the method of the calculations is primarily determined by the ratio h/ R. This can be understood 
by looking at the single-particle unperturbed potential (see Fig. 2). 
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FIG. 2. The cross section y = Q of the single-particle potential V\(x,y). 
and the solid line is the Hydrogen-like molecule {hi, /12 < R). 



The dotted line is the Helium-like atom {hi, /12 > R), 



Function Vi{x, y) has either one or two minima. The transition from one case to another happens when a minimum 
near the origin changes for a maximum. For the case of equal heights hi = h2 — h the function Vi{x,y) has one 
minimum if h/ R > l/\/2- 

The case ft.i, ft.2 < R corresponds to a 2D Hydrogen-like molecule, when there are two well-separated minima. In 
this case the splitting is exponentially small and determined by the degree of overlapping of wave-functions in the 
classically forbidden under-barrier region (especially at ri,r2 ~ 0) . 

The opposite case /ii,/i2 > R corresponds to a 2D Helium-like atom when only one minimum exists, and the 
interaction term l/ri2 should be taken into account from the very beginning. 

If the difference between hi and /12 is large enough then the electrons will spend most of their time in the potential 
well of the impurity closest to the plane. The potential Vi{x,y) has one minimum in this case. Therefore this case 
corresponds to a perturbed "Helium" atom (or perturbed ion H^). 

The paper is organized as follows. In Section 2 we consider the case of a hydrogen-like molecule {R > 2h). 
By making use of semiclassical methods we obtain analytical results for a pure two-dimensional hydrogen molecule 
(/ii, /12 = 0, i? ^ 1) 1^. The perturbation theory results (/11.2 <^ I < R) and the saddle point methodfl < ft.1.2 ^ R) 
represent other analytical limits in this case. The details of the calculations are given in appendices |^,^ 

In Section 3 we discuss a helium-like atom (/ii,2 > R)- We start from the result for pure two-dimensional helium 
PI (i? = /11.2 = 0) and calculate the perturbation theory corrections. The quadratic approximation is applied to the 
limit hi^2 > R > I- Section 4 is devoted to the interpolated formula, which connects the two different limits for the 
distance R at fixed vertical coordinates of impurities /ii, /12. Our results are summarized in the Conclusion. 



II. HYDROGEN-LIKE MOLECULE 



A. General formula, R^ 1. 



The exchange constant Jik is defined as one-half of the energy difference between the lowest two-particle energies 
for total spins S = 1 and S — 0: 
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Here ^1^5 and '9 a denote the lowest symmetrical and antisymmetrical two-particle coordinate wave functions of the 
singlet and triplet states respectively. When /ii 2 < R and i? ^ 1 the most appropriate method to determine the 
term splitting due to the spin-spin interaction is the Gor'kov-Pitaevskii method [Q-|6|. We apply it here to the two- 
dimensional electron system. In accordance with this approach we reduce the expression for J to a surface integral 
over a hyperplane in the 4-dimensional coordinate space. To this end, we multiply Eq. (Q) by ^^ and Eq. (^) by ^5, 
take the difference between the results and calculate the integral over a part of the volume in the four dimensional 
configuration space of the electrons. We choose the integration volume as a region xi < X2- This is a four-dimensional 
volume to the left of the hyperplane E determined by the equation xi — X2. Using the Hamiltonian (|l|) and Gauss' 
theorem we obtain 
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where cr is a closed hypersurface, and Q is the volume bounded by a. 

Now we introduce a combinations of the functions ^1.2 = l/\/2(^s ± '^a)- If the phases of '^s,A are properly 
chosen, the function ^1 will be large only when electron 1 is near the left minimum and electron 2 is near the right 
minimum. The function ^2j obtained by permuting ri with r2, is localized almost entirely on the "far" side. A simple 
calculation gives 



*s^A dri dra = - / / (*i 



*^)dridr2~ 1/2. 



Substituting the wave functions ^1.2 into Eq. 
over the hyperplane S. 



2J = -2 



dxi 



we finally obtain that {Es — Ea) can be expressed as the integral 
dx2 dyi dy2- (7) 
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dxi 



Eq. (|^) shows that the main contribution to the exchange constant is given by the region where the electrons are 
close to each other. Indeed, the x coordinates of both electrons coincide (xi = X2), however, the y coordinates may be 
different. In this case there are strong correlations between the positions of the electrons due to Coulomb repulsion. 
This means that we should go beyond the approximation where the two-particle wave function of the electrons is 
represented as a product of single-particle wave functions. 

The details of the calculation of the wave functions ^1.2 and the surface integral (M) are presented in the appendices 
0,g. The final result is: 



2 J(a, /3, R) = R-+Ti--e-^"^ [D {a, I3,R) + D (/?, a, R)] , 
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where the function D {a, /?, R) is determined by Eq. ( B12| ). Here 
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are "ionization energies" in the potential I/a/t^+T?, and (i = a + (3. 



B. Two-dimensional hydrogen molecule and perturbation theory for hi, /i2 «C 1 

In the case a — (3 the function D{a, a, R) is independent of R: 

Do{a) ^2D{a, a, R) = SV^A^ (J-^ T (^^^ ^' exp {-t/a) (2 - t^/^'^t^/'- dt. 



Therefore for the case hi = /12 



2J(a, a, R) = Do{a)R'/'^°' exp(-2ai?) 



(9) 
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Here Ac is the normalization constant for the single particle wave function in the potential 1/ \/r^+h? (see Appendix 

The case h — Q corresponds to the two-dimensional hydrogen molecule. In accordance with Appendix [Aj we 
substitute a = 2 and A^ = A/\/2t: in Eq. (|lO|) and we obtain 

2 J(2, 2, R) = 30.413 R^/^ exp(-4i?) (11) 

We can carry out further analytical estimates for small ft.i,2- When h is much smaller than the Bohr radius we take 
Aa from (A£), put it into Eq. (0) and obtain 



D^{a) - ^^ "-' ' ^ \ exp (~t/a) (2 - t)^/^" t^/s. ^^ 
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Function I?o(q^) is a regular function near the point a = 2 (see Fig. 3). It can be expanded in a Taylor series: 

Di^{a) « 30.413 (l + 1.665(a - 2) + 0.456(q - 2)^ + . . .) (13) 
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FIG. 3. Exact (crosses), linear (dashed line) and quadratic (solid line) approximations for function Do (a) 
Taking into consideration that (see Appendix H) 

a w2- Aij/2 
with 

Ab = 16/i(l + 2/iln(2/i)), 
finally we obtain 



J(a, a, R) « J(2, 2, R) [1 - 0.833Ab] R^^'' expiAsR) w J(2, 2, R) 
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Thus, the sphtting increases with h. Indeed, when the impurity centers go out of the plane the interaction between 
them and the electrons becomes smaller. Therefore, radii of bound single electron wave functions and overlapping 
increase. 

When hi 7^ /12 and /ii, /i2 ^ 1 the formula for J should be rewritten in the form: 

J{a, /3, R) « J(2, 2, R) [1 + Asi-di +R + 7/16) + (J^ds] , (15) 

where 

Ab,,5b = (Ab, ±AbJ/2 

and di, ^2 are some coefficients which can be extracted from the expansion D{a, (3, R) for small /ii,2- 

C. 1 < /i < i?. Saddle point method. 

At /i ^ 1, one can easily find that E{h) w —1/h (See also Appendix ^). Then a w \/2/^/h <C 1. and the integral 
in -Do(ct) Bs well as the Gamma- function can be estimated by the saddle-point method: 

r(A/2 + l/2)«V^KA^(^) ' °^P(-^)' 
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where A = 1/a. Thus, 
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and 



2J(a,a,i?) = W — -iU^Se^l -'/('") i?7/(2a)exp(-2ai?), /i » 1, i?//i » 1. (18) 



III. HELIUM-LIKE ATOM 
A. Pure 2D Helium atom. 

When two impurities are close to each other (R ^ \/2h) there is only one minimum for the single-particle potential 
Vi{x,y). Hence, the semiclassical expression ([7|) for the exchange constant is not applicable anymore and we have to 
exploit other methods. 

For the case /11.2 = 0, R ~ we use a variational approach for the two-dimensional helium atom. We consider 
impurities as the "nucleus" with a charge Z — 2 and the trial functions are 

*A = e^°" [at coshipt) + {as - 2) sinh(/3i)] , 

s,t = ri±r2. (19) 

Here a and /3 are free parameters, which have to be determined from the variational principle. The details of the 
calculations were presented in ||^. We obtained the following results for the energies: 

ECS) = -11.635 a.u. 

E{^S) = -8.193 a.u.. (20) 

Thus, in this case 

2J = 3.567 (±l%)a.u. (21) 



B. Small corrections due to finite size of the "nucleus" 



When the distance between impurities (or distance from impurities to the plane) is smaU in comparison with the 
Bohr radius we can use first-order perturbation theory to estimate the level's energy shift: 






But the integral over all variables except one is equal to 
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where p(r) is the electron density (normalized to two - the number of particles). Therefore, 
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(22) 
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Here we take into consideration the fact that the main contribution from density gives the term which is simply 
proportional to the exponent (for the S = state ao — '^Zq and for the S* = 1 state ai = Zi + Z2/3). 

For the case hi ^ h2 = 0, R ^ 1 the difference between the point "Helium" nucleus and the real potentials for the 
single electron is 



AF(r) = - - 
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r yjr^ + a^ — 2ra cos(0) \J r^ + a^ + 2ra cos(0) 
The double integral for the last two terms of (E3) can be reduced to the single integral 
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Here K{x) is a complete elliptic integral of the first kind. 

For an estimate of the first integral we can use the fact that exp(— 2aQ;x) ~ 1 in the region r < a and in the second 
integral we put 
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Thus 

A^ = 8ap(0)(l + ln(2)), 
and hence, the exchange constant is 

2 J(0, 0, R) = 2 J(0, 0, 0) - 4 (1 + ln(2)) Ap{0)R = 3.567 - SOAR. 
Here 



Ap(0) = po(0)-pi(0) 



For the case R — and hi , /12 <C 1 
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Therefore, 



J(/ii, /i2, 0) ~ 2 J(0, 0, 0) - 27r {hi + /12) Ap(0) = 3.567 - 28.2(/ii + ^12). 



(30) 
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It is worthwhile to stress that the convergence region of the perturbation theory is very small (i? < O.las) because of 
the large numerical factors in the second terms of Eqs. (ESh and (p^) . 



C. Case R <C h^'"^, ft ^ 1. Quadratic approximation 

If the impurities are far enough from the plane hi ^ h2 — h ^ 1 and \/2h > R, the single-particle potential (0) has 
an oscillatory shape near its single minimum. Therefore, taking into consideration its Taylor expansion near zero the 
potential in two-electron hamiltonian (Q) has the following form: 

u{xi,x2,yi,y2) = -^Tp^p + ^,^2 + Jf^ [{^l + 4 + yl + yl] + 



{h^ + a^f'"' "' ^{xi - x^Y + (2/1 - y2) 
Introducing the notations: 
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and substituting the variables 
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we separate center- mass (^, m) and relative motion (77, u): 
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The energy splitting between the singlet and the triplet states is determined by the relative motion part of the 
Hamiltonian. If we neglect by the last "perturbed" term in (p4h, the potential for the relative motion is independent 
of Lp. Therefore 



^(77, w) = (/)(r) exp(zr7i(^)/v 27r. 

The state with 5* = corresponds to an even function of coordinates and the state with 5 = 1 corresponds to an 
odd function. Parity is determines by the two-dimensional angular momentum m. 

r —t — r => r —I- r, (^ — J tt + (/?. 

Thus, symmetric (S* — 0) wave functions correspond to even angular momenta (|7ti| ~ 0,2,4...) and antisymmetric 
wave functions correspond to odd angular momenta (|7tt-| = 1,3, 5...) 

The exchange constant is 2J = £^01 ~ -E'OOj where £^„^jm| are the energy levels for the potential 

U{r) = 4r'+ ' 
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This potential can be approximated by an oscillatory one after expanding around its equilibrium position rg 

2-1/6/4/^: 
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Thus, in this limit we obtain 
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2J«V3c.+ =(^9i^^-^^^j , i?«/^. (35) 

Taking into account that first-order corrections must be small in comparison with the energy splitting we obtain more 



accurate condition for the applicability of Eq. (35) 



< ^fr ^ - ^ « 1. (36) 

It is equivalent to the condition R <^ h?/^. 

IV. INTERPOLATION 

Let us consider equal heights of the impurities above the plane: hi = h2 = h. In the asymptotic region, for large 
distances between the impurities, we have 

2J{h, R) = L>o(a)i?^/2" exp(-2ai?), R:>h. (37) 

Here Do{a) is determined by Eq. (||) and — [a(^)] /2 is the single-particle ground state energy in the potential 
— (r^ + /i^)^^/^ (see Appendix ^). In the opposite helium-like limit we have 

r -3.567, h = 
2J{h,0) = \ AT ■ <, (38) 

Let us now make a plausible assumption that for small distances between the impurities {R ^ h) the behavior of the 
exchange constant is 

2 J(/i, R) « 2 J(/i, 0) exp {-jR - 72 i?^) , (39) 

where 7 = j{h) is a free parameter, and 72 = 72(7)- Such form of the functional dependence follows from the fact 
that the most accurate numerical calculations for the three-dimensional hydrogen molecule H are fitted accurately 
at small distances by the formula 2 J(0) exp(— 7ii? — 72 -R^) with just two independent parameters 71, 72 g. Let us 
match the second derivative of In (2 J(/i, R)). In the two regions it must have the following behavior [See Eqs. ( p7tp9| )] 

52ln(2J) f -272, R<h 

^' 7_ 

The simplest formula that satisfies both conditions is 

dHnjJ) _ 272 

(9i?2 ^ 1 + 4/772ai?2 

After integrating over R twice we obtain 

ln(2 J) = C -jR AR arctan(Ai?) + — ln(l + A^R^). (42) 

2a 4a 



dR^ "1 -^7^, i?»/i. ^^°^ 



(41) 



where C and 7 are integration constants and A — ■y/4a72/7. It is obvious that we have to put C = In {2J{h, 0)) in 
order to fulfill (p3) . On the other hand to guarantee the correct asymptotic exponent of Eq. ( pT\ ) we have to impose 
a constrain 



^^4a(2a-7)^ (43) 
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Finally, the interpolated formula has the form 
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2J{h,R) = 2J{h,0)[l + A^R'^]' exp 
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(44) 



Two parameters in this formula, J{h, 0) and 7(/i), are not defined yet. Additional numerical calculations are needed 

For 
, to the 



to determine them. The parameter 2J{h,0) decreases from 3.567 at ft. = to \f%jW for /i ^ 1 [see Eq. 
intermediate values of h it can be determined by applying a variational approach similar to that of Ref. 
helium-like atom. We can do a rough estimate of J{h,0) using an interpolation between these two limits: 



2J(/i,0) 



3.567 
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The parameter 7(/i) can be estimated by considering small and large h. In our previous paper [^ we have found a 
reliable interpolated formula for splitting in the case h = 0. It reads 



2 J(0, R) = 3.567 (1 + IMR^y^^ exp {-0.3R - 2.355i?arctan (1.346i?)) . 



(46) 



To get this equation from Eqs. (0),(^) at ft ^ one should put 7(0) — 0.3. We can argue that 7(ft) should decrease 
with increasing ft. Indeed, Eq. ( |35| ) shows that function J(ft, R) is an algebraic rather than exponential function of 
R at R <$:. h^/'^. To satisfy this condition one should assume th at 7( ft) ^ h^^/^ at ft ^ 1. Then it follows from Eq. 
(H) that A = %c? jl-i: at ft > 1, where a is determined by Eq. ( |A10| ). Thus A ^ l/h at ft > 1. 
The function 7(ft) can be chosen, for example, in the following form 



7(ft) = 



0.3 



1 + ft3/4 ■ 



It takes into consideration the main features. We would like to stress that there is practically no dependence on 7 in 
the exchange constant (H) , since it is small even for ft = and it disappears completely due to cancellation of two 
terms in Eq. (|4J) at large R. 

The results are summarized in Fig. 4, where a semi-logarithmic plot of the exchange constant vs. R for three 
different heights is presented. The solid line shows 2 J(0, R) for the pure two-dimensional hydrogen molecule. The 
curve ft = 0.8 corresponds to the atomic parameter a « 1 and, consequently, it has a two times smaller asymptotic 
exponent in the exchange constant. Of course, this case corresponds to the intermediate values of the parameter ft. 
Therefore, an approximate realistic magnitude of 2 J(ft, 0) was chosen. The last curve shows that in the limit of large 
ft the behavior of the splitting energy is flat, and almo st independent of R for R <C ft^^'*, in accordance with Eq. (|3^). 
The value of a is taken from the oscillator limit ( |A1G| ) and 2J(ft,0) corresponds to the lower formula in (|3q). 
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FIG. 4. The exchange constant as a function of R and h. The curves for three different heights are presented. The crosses 
correspond to the numerical data for _R = and R — 1 and asymptotic form ( |ll| ) for the pure two-dimensional case {h — Q). 
The solid line is the interpolated formula (k6|). Dot-dashed and dotted lines present the interpolated formula (l44) with h = 0.8 
and 10 respectively. 



V. CONCLUSION 

We have calculated analytically different limits of the singlet-triplet energy splitting in the two-electron Hamiltonian 
(|l|) in the following cases. 

1. Hydrogen-like molecule. 

(a) /ii, /i2 = 0, i? > 1, which corresponds to a 2D Hydrogen molecule, Eq. (jll|). 

(b) /ii,2 <C 1, i? > 1, a Hydrogen-hke molecule, where perturbation theory is applicable, Eqs. (||,|o||lJ,|l^). 

(c) 1 <t: hi,h2 <t: R, saddle point method, (Eq. (|l|)). 

2. Helium-like atom. 

(a) R = hi = h2 = corresponds to the 2D Helium atom, Eq. (pl|). 

(b) i? <C 1, hi ^ h2 — 0, 2D Helium atom with perturbation theory corrections, Eq. (|28|). 

(c) R — 0, hi, ft,2 ^ 1, 2D Helium atom with perturbation theory corrections, Eq. (p^). 

(d) hi,h2 > R > 1, quadratic approximation for the Helium-like atom, Eq. (p5|). 

Finally, we have presented the interpolated formula (Q) connecting the limits 1 and 2. This formula can be used 
for calculations of the ground state in spin glasses and semiconductors based on the Heisenberg Hamiltonian. We 
have shown that the i?-dependence of the splitting becomes weaker with increasing /ii, /i2. 
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APPENDIX A: EIGENVALUES AND EIGENFUNCTIONS IN 2D CENTRAL POTENTIALS 

The Schrodinger equation for a radial wave function in a central field is 



9^ Id m? 

Qj.2 J, Qj, J,: 



l^,+Z—- — -'^U{r) + 2E\R{r)^Q. (Al) 



For the Coulomb potential, U{r) = —Z/r, the eigenvalues and eigenfunctions are well-known: 

I"! 

Er, = -, TTTT, (A2) 

2 (n-- 1/2)2' V ; 

n= 1,2,3... 
TO = — n + 1, —n + 2, ...n — 1 

Ze*™"^ / 16 (n- |to| - 1)! r '^-^ "1 '"' 



•^rvmir) = 



2tt y (2n- l)3[(n+|TO| - 1)!] 



4rZ 
2n-l 



exp ( -Ji:^) L^l™l , ^ (^^ 1 , (A3) 



where Lf{p) is the generalized Laguerre polynomials, and the level degeneracy is g{n) = 2n — 1. 
The most important for applications here Is and 2s states have the following form. 

\/ Zn 
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Thus, for the ground state with Z = +1 
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Let us consider now the potential 



C/(r) = -1 



/l2. 



(A6) 



Far away from the "atomic residue" U(r) ~ — 1/r and the wave function for the s-state obeys the Schrodinger equation 

A 1 o? 



with the solution 



V' 



(r) = A^ri/^-i/^e""'- 



(A7) 



up to (/j/r^ accuracy. Here Aa and a are atomic parameters. Their magnitudes are determined by the behavior of 
the electron inside the "atom" . 

The results of the numerical calculations for the normalization constant and the eigenvalues are presented in 
Figure 5. 
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FIG. 5. (a) Ground and first energies as functions of the parameter h. Solid lines are numerical results. The dot-dashed 
lines correspond to perturbation theory for small h and an asymptotic oscillatory energy ~h~^ + h~^'^ for /i 3> 1. (b) The 
constant Aa for the ground state. Numerical results (solid line) and analytical es timat es for both regimes are shown. Dotted 
line corresponds to Eq. (A.9) and dot-dashed line is the quadratic approximation (A.12). The inset is a log-log scale plot. 



Simple analytical estimates for the cases of small and large h for ground state energy and the asymptotic coefficient 
Aa can be made. 

When /i <C 1 the potential is only slightly different from the Coulomb potential and the contribution to the energy 
can be obtained using perturbation theory: 



A£-(0|V^|0) = 16y°°('i 



1 



V/i^ + r"^ 



exp(-4r)r dr « 16/i (1 + 2h ln(2/i)) 



(A8) 



Note, that due to the big factor in front of fe in (A8) the convergence radius for the perturbation series is small enough. 
Supposing that the asymptotic form (A7) of the wave function is valid for all values of r, we obtain the following 
simple estimate for Aa 



2-KAi I r^/" exp(-2ar) dr = 1 ^ A^ 
'o 



2l/a-l/2^1/a+l 

V'^r(2/a) 



(A9) 
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In the opposite case, when h^ 1 the solution has to be close to the oscillatory one 

1 1 



E = -- 



h /i3/2 



(^ + 1), 



(AlO) 



where n = 0, 1, 2.... The main contribution to the normalization comes from the Gaussian part of the wave function 
because the contribution of the Coulomb tails is negligible. One gets for the state with n — 



where uj — l/Zi^/^. Then, 



ip{r) w iposc{r) 



Aa — 



y^ exp (-wrV2) , 



0F/i3/4 ■ 



(All) 



(A12) 



The agreement between the behavior of the estimate (A12) and the numerical calculations can be observed in the 
inset in Fig. 3. 



APPENDIX B: EXCHANGE CONSTANT FOR 2D HYDROGEN-LIKE MOLECULE 

Since J{R) is exponentially small as i? ^ cxd, ^a and ^s are solutions of the same Schrodinger equation, and 
therefore, with exponential accuracy their combinations 



*i,: 



V2 



are also the solutions of the same Schrodinger equation with the Hamiltonian (|l|). They correspond to the states of 
"distinguishable" particles, when, e.g. for \['i(ri,r2), the first electron is principally located near first ion at x = —a 
and the second electron near the second ion with x = a. The electron energy 

is accurate to terms '^ 1/R^. Therefore we are looking for \l'i^2 in a form 

'^i{ri,f2) = (j)a{\ri + a\)(j)^{\f2 - a\)xiri,r2), 

*2(ri, fa) = (/)„(|f2 + a\)(j)^{\fi - a|)x(f2, fi), (Bl) 

where (j)a,i3 have an asymptotic behavior (see Appendix |A|) 

(^„(r)=A„rl/"-l/2e-"^ 

and X is a slowly varying function of rl, r2. Substituting $i into wave equation and neglecting the second derivatives 
of X, we obtain 



^ dx g dx 
dxi dx2 



1 



1 



1 



1 



^{xi -~X2Y + {yi - J/2)2 2a a -XI a + X2 



X-0. 



(B2) 



Equation (B2) valid under conditions 



\xia\ < a, 2/12 = \yi ~2/2| < Va 
Ra,RI3-» l,R\a- (3\ < 1. 



The general solution of (B2) is 

F{Ci(xi,X2),C2{x,Xl,X2,Vl2) = 0, 

where Ci, C2 are integrals of the motion of the ordinary differential equations: 



(B3) 
(B4) 



dx\ 

a 



dX2 



dx 

X 



1 



1 



Vi^i -XiY +2/12 ^" a-xx a + X2 
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Hence 



dxi dx2 



a 

dxi 
a 



f3 



Ci = xi/a + X2//3, 



C2^X 



y^Jxr^^X2f^^2 ^" a-xi a + X2 
exp(2:i/2aa)[a - xi]i/"[a + X2Y''^ 



= -d\n{x) 



y/{xi - X2Y +2/12 - Xi + X2 



Combining (B4p3E|) we can write the function x a-s 



X{x\,x2,y\2) 



exp(-5l^) \J{X\ -X2Y^vl2-Xx +X2 



[a-xi]i/"[a + a;2]i//3 



-/ 



X2_ X2_ 

a f3 



(B5) 



(B6) 



(B7) 



where unknown function /(u) is determined from the fact that x — ^ 1 when xi — > —a, X2 is arbitrary, or when 
X2 — > a and xi is arbitrary. Finally, after expanding |r ± a| ~ |a ± a;| + t/f2/2|oi a;| in the exponent, we obtain 



2-13 



'i'iiri,f2)=AaAi3{a + xi) 2° (a-a;2)^^exp 
X{xi,x2,yi2) = 



-a{a + f3) + (3x2 — oixi 



ayl 



Pvl 



2{a + xi) 2(a — ^2) 



X(ri,f2), (B8) 



e 2a 



2a 



a— xi 
2a 



a+a:2 






/3" 



/3(a+a^i)+a(a+a;2) 
a+X2 

/3(a— a:i)+a(a— 3:2) 



\J (x\-X2Y-\-y\^-x\-^X2 



^/(/3(a+a;i)+Q(a+a:2))^ + (Qyi2)^+/3(a+3:i)+a(a+a;2) 



■y(xi-X2)2 +1/^3-3:1 +2^2 



V(/5(a-^l)+Q(a-a;2))H(/3yi2)2+/3(a-xi)+Q(a-a:2) 



, (B9) 



» + /3 



Here the upper expression is given for x\+ X2< 0, and lower expression for xi + 2:2 > correspondingly. 
Substituting (B8) in Eq. (H), and differentiating only the exponent we obtain 



2J = +2(a + /3) / [*i*2],,=,,^, dxdyidy2 

J —a 

Introducing notations fi = a + f3 and u — f3 — a and taking into consideration that in the approximation (B3) 



(BIO) 



yj[{0 + a){a-x)f + {(3yi2y + (/3 + a){a - x) ^ 2(/3 + a)(a - x). 



the formula ( |B10D finally transforms to the 



2 J(a, P, R) = Ri+^-T^e-"'^ [D (a, p,R) + D (/3, a, R)] 



where D (a, /3, R) is the following function 



.,o,,,«,=4^.^4(f)-"ri±^ (2-..'.3- 



i 



2/i y V a 

1 exp (-(1 - a;)/a - vRx) (1 + 2;)2//3-2/a+l/p(^ _ .^)2/a-l/A< 



(l-.Tl./Ai)'+'/'' 



dx. 



In the case a = /3 it is independent of R at all 

L>o(a) = 2i:'(a,a,i?) = 8V^A;^ ( ^ ) T 
In this case 



l\'/'"„/a + l^ '^ 



4a / V 2a 



exp 



(-t/a)(2-t)l/2a^3/2a^^_ 



(Bll) 



(B12) 



(B13) 
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2J(a, a, R) = Dfi{a)R^''^°' exp(-2ai?) (B14) 

For 2D H2 molecule {a ==2, Aa ^ 4/\/27r) it gives 

2 J(2, 2, R) = 30.413 R'/'^ exp(-4i?) (B15) 

For comparison we remind that in 3£'-case 

2J3D = 1.636i?^/2gxp(-2i?) 
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